Unified quantitative model for magnetic and electronic spectra of the undoped 

cuprates 
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Using low-energy projection of the one-band f-t'-t"-Hubbard model we derive an effective spin- 
Hamiltonian and its spin- wave expansion to order 1/S. We fit the spin- wave dispersion of several 
parent compounds to the high-temperature superconducting cuprates: La2Cu04, Sr2Cu02Cl2 and 
Bi2Sr2YCu208. Our accurate quantitative determination of the one-band Hubbard model param- 
eters allows prediction and comparison to experimental results of measurable quantities such as 
staggered moment, double occupancy density, spin- wave velocity and bimagnon excitation spec- 
trum and density of states, which is discussed in relation to K-edge RIXS and Raman experiments. 

PACS numbers: 74.72.Cj, 75.30.Ds, 78.70.Ck 



High- Tc superconductors challenge all known theoreti- 
cal approaches by mixing charge and magnetic degrees of 
freedom and lacking a small variational parameter. The 
One-Band Hubbard Model (IbHub), proposed by Ander- 
son [l| to describe their Cu02 planes, includes both of 
these difficulties in its parameters, the electron filling, the 
hopping matrix element t and the Coulomb repulsion U . 
The ratio t/U is moderately small in the cuprates, thus 
most approaches start by projecting out double occupan- 
cies (DO) to obtain the Heisenberg model at half-filling 
or the t-J model for hole- or electron-doping. Such a pro- 
jection is too complicated to be carried out exactly but 
may be performed as an expansion in powers oi t/U 

Experimentally, different techniques have probed sepa- 
rate channels - magnetic or electronic. Inelastic Neutron 
Scattering (INS) on La2Cu04 (3 demonstrated that the 
projection must be carried out at least to the fourth order 
(t^/U^) to reproduce the observed magnetic excitation 
spectrum. Angle-Resolved Photo-Emission Spectroscopy 
(ARPES) indicates that first, second and third nearest 
neighbor hopping matrix elements are needed to repro- 
duce the observed electronic quasiparticle dispersion Q. 
A quantitative description of the undoped high- Tc parent 
compounds therefore needs a model incorporating both 
those considerations 

In this letter, we develop such a quantitative the- 
ory and present the resulting sets of IbHub parameters 
for single and bilayer undoped cuprates Sr2Cu02Cl2 0, 
La2Cu04 and Bi2Sr2YCu208 [this work], obtained 
by fitting their magnetic excitation spectra measured by 
Resonant Inelastic X-ray Scattering (RIXS) and INS. 
We start from the IbHub Hamiltonian given by 
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FIG. 1. (Color Online) Cu02 planes from the perovkite struc- 
ture. Examples of exchange loops from the effective spin 
Hamiltonian of Eq.[2]are indicated with first, second and third 
nearest neighbor hoppings t, t! and t" . 



sites i and j, \J the effective on-site repulsion and ni^a = 
c\ „Ci^a the fermionic number operator. At half filling, the 
kinetic term mixes states with different number of DO, 
which in the limit t/U <^ 1 separate into different energy 
scales. We use the unitary transformation technique ^ 
up to order t'^/U^ to decouple these states and obtain 
the effective spin Hamiltonian 
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where i?*^^^ is a constant and {f^}, |*^| ^^^^ 
stand for the plaquette ensembles of two, three and four 
sites respectively, and connected as sketched. We em- 
phasize that the loop ensembles are fully defined by the 
where cj ^ and Ci^cr creates or destroys a fermion with spin considered hopping matrix elements iy- and the lattice 
a on site i, tij is the hopping matrix element between topology, so that the Hamiltonian of Eq. [5] is the gen- 



2 



eral low-cncrgy projection to order t'^ /U^ of the IbHub 
at half fining for any lattice. Examples of plaquettes of 
two, three and four sites involving first, second and third 
Nearest Neighbor (NN) hopping amplitudes t, t' and t" 
are sketched on Fig. [TJ A similar development but away 
from half- filling would result in a t-J model with the 
same magnetic couplings as in Eq. [2] plus a family of 
charge plaquette hoppings. Being the canonical model 
for high-Tc superconductors, a large variety of analytic 
and numerical approaches exist to study such a model 
i. 

In contrast to the doped compounds, the derivation of 
measurable quantities is much easier for their undoped 
parents. Here, we use Spin- Wave Theory (SWT) to de- 
rive the dispersion of the magnetic excitations which arc 
measured by INS or RIXS. We expand the spin operators 
of Eq. [2] in terms of Holstein-Primakov bosons. Keeping 
the first 1/5' correction for the tfj/U terms, the Hamilto- 

nian transforms as TL^^^ = EN+n2+'H4 + 0{l/ S), where 
the Neel ground-state energy E-^ includes the constant 
of Eq. [21 and the harmonic dispersion luq (k) is obtained 
from a Bogoliubov transformation of the quadratic term 
7^2- A Hartree-Fock decoupling Q of the quartic term 
T-Li results in an overall momentum-dependent correction 
to the magnon energy so that our final spin- wave Hamil- 
tonian reads 

H^*' = ^c(k)LJo(k)4ak + En+5E (3) 

k 

where a's are free magnon operators, Zc{^) is the 1/S 
renormalization of their dispersion and SE is the quan- 
tum correction to the ground state energy at this order. 
For the bilayer square- lattice, Eq. [2] is still valid but the 
loop ensembles now include inter layer hopping t±, and 
two boson flavors to account for the top and bottom sites. 
This results in two magnon modes which are gapped re- 
spectively at (0,0) and (7r,7r) but degenerate along the 
magnetic Zone Boundary (ZB). 

Before applying the result to the cuprate compounds 
SCOC and LCO 0], we report here new measurements 
of the spin-wave (SW) dispersion in Bi2Sr2YCu208, a 
bilayer parent compound. Single crystals were grown by 
the flux method with yttrium ensuring an insulating an- 
tiferromagnetic phase. The SW dispersion was measured 
using Cu edge RIXS at the SAXES end-station of 
the Swiss Light Source ADRESS bcamlinc, experimental 
details and data analysis as described previously 

The spin-wave dispersions of the various compounds 
are shown in Fig. [2l They all feature a dispersion be- 
tween the ZB points (tt, 0) and (7r/2,7r/2) which can in 
principle be explained by the effective model of Eq. [2] with 
NN hopping alone (dashed red lines). However, this ap- 
proach results in unphysically low U = 2.2 eV 0] for LCO 
and J7 < 2 eV for SCOC and BSYCO. From the former 
to the latter, their ZB dispersions respectively reach 40, 
70 and 55 meV. Although U is an effective on-site repul- 
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FIG. 2. (Color Online) (a) Fits of the measured SW disper- 
sions (solid lines) and the IJHei SW dispersion (dashed lines) 
with Jnn = a;(.^/2 7r/2)/2^c. From top to bottom: BSCCO 
(Jnn = 0.15 eV), SCOC (Jnn = 0.12 eV) and LCO 0] 
(Jnn = 0.14 eV). (b-d) Same as in (a) but subtracted from the 
IJHei SW dispersion to enhance details, BSCCO, SCOC and 
LCO respectively. Data points from experiments are folded 
onto the equivalent high-symmetry axes of the first Brillouin 
zone. 



sion, closer to the charge-transfer gap than to the bare 
Coulomb repulsion, a good IbHub must use an effective 
parameter compatible with electronic and optical spec- 
troscopies, which request U ~3-4 eV for the cuprates [lO|. 
Consistently with ARPES results, we therefore include 
second and third NN hopping in our effective model and 
derive the spin- wave dispersion of Eq. |3] which is now a 
hmction of four parameters {U,t,t' ,t"). The measured 
SW dispersions contain three distinct constraints, the 
(7r,0) and (7r/2,7r/2) ZB energies and the spin-wave ve- 
locity. We thus expect a one-dimensional solution and 
choose the free parameter to be the effective U. The fit- 
ting procedure is as follows. For a fixed choice of {U, t" /t) 
we start by fitting the two other parameters t/U and t' /t. 
As the calculation of 1/S estimate of Zc involve a slowly 
convergent integration over k-space, we include it in a 
two-step iterative approach. First, we fix its value to 
the uniform Zq = 1.1579 obtained for the NN Heisen- 
berg model (IJHei). Then, we fit ZoUJk{U,t/U,t' /t,t" /t) 
using a non-linear least-squares algorithm and calculate 
a first non- uniform Zi (k) from the obtained parameters 
set. We iterate this procedure until Zn{^ converges, typ- 
ically after 10-15 steps. In the case of BSYCO, we further 
include an interplanc hopping However, the resolu- 
tion of RIXS does not allow to distinguish the splitting 
between the two magnon modes and we fix it to the value 
t j_ = 54 meV reported by Chuang ei al. [llj . 
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FIG. 3. (Color Online) Summary of the fitting results, (a-c) Goodness (x^) of the {t/U,t' /t) fit for fixed (U,t" /t) in the case 
Bi2Sr2CaCu2 08, Sr2Cu02Cl2 and La2Cu04 respectively. Solid (blue) and dashed (green) lines are best fit solutions as function 
of U for t't" < and t't" > respectively, (d-f) along the best fit lines defined in (a-c). (g-i) {U,t,t' ,t") solutions along 
the best fit lines defined in (a-c). (j-1) Staggered magnetization, DO density and k-averaged quantum 1/S renormalization Zc 
along the t't" < best fit lines defined in (a-d). Lines: Solid blue BSYCO, dashed red SCOC and dotted black LCO. Shaded 
areas indicate k-variation of Zc 



The fitting results over the {U,t"/t) plane are shown 
in Fig.Owitli the {U, t/U, t' /t, t" /t) parameters along the 
best fit lines. Overall, the four compounds share common 
features i.e. a strong lower boundary for J7, an increase 
of \t'\/t and t" /t with U , and a slowly varying t/U . Due 
to the t^t't" /U^ term, the calculated magnon dispersion 
is symmetric in the signs of t and t' but not in the rela- 
tive sign of t' and t" resulting in two separate solutions 
for t't" < and t't" > 0. From the best-fit lines, one 
can see that the inclusion of \t"\ is necessary in order to 
get U « 3-4 eV. For some regions of the {U, t" /t) space, 
the Neel state is not the classical ground state of Eq. [51 
and/or it is destroyed by quantum fluctuations. Both 
cases can be systematically determined by looking at the 
size of zero-point fluctuations {a\ai). The outermost void 
regions in Fig. [3Ia)-(d) are those where Neel order is un- 
stable. In Fig. E^m), we calculate the evolution of the 
staggered magnetization as a function of U along the 
best fit lines. Increasing [/, t' and t" grow while t/U 
stays roughly constant, bringing more frustration and 
subsequently reducing the ordered moment. We calcu- 
late the double occupation density using the Feynman- 
Hellmann theorem {rii^^ni^]^) = d{w'^^) /dU . Along the 
best fit lines, a ?7- independent value 5% is found for all 
cuprates in agreement with the electronic shielding fac- 
tor calculated in ref. [l2|. The l/S estimate of Zci}f.) is 
found to vary only about 2% across the Brillouin zone. In 
Fig. El^o) we show the average value, which vary from 1.2 
to 1.3 for U K, 3-4 eV which again reveals the prominent 
role of quantum fluctuations in the range of parameters 
relevant for the cuprates. 

The effective U cannot be directly obtained through 
magnetic excitations. However, more direct cxperimcn- 
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BSYCO 


470(10) 


-205(3) 


79(4) 


0.3 


0.146 


5.9% 


SCOC 


480(10) 


-200(5) 


75(5) 


0.29 


0.163 


5.1% 


LCO 


492(7) 


-207(3) 


45(2) 


0.4 


0.195 


5.2% 



TABLE I. Parameters determined from the spin-wave fits 
with fixed (7 = 3.5 eV and t't" <0. 



tal techniques may give good estimates of U thus de- 
termining a unique set of IbHub parameters for each 
of the above compounds. In particular, an estimate 
of the IbHub parameters of SCOC was obtained from 
ARPES as U = 3.5 eV, t = 0.35 eV, t' = -0.12 eV and 
t" = 0.08 eV. Consistently with those parameters and in 
order to compare the three cuprate compounds, we adopt 
in Tab.Ua uniform value U = 3.5 eV and the t't" < so- 
lution. A more accurate determination of U could be 
found in the charge-transfer (CT) excitation part of the 
RIXS spectrum. Using the above ARPES parameter esti- 
mates, Hasan et al. could identify a dispersing excitation 
around 3 cV in Cu K-edge RIXS as CT excitation [l3l |. 
A similar approach using our parameter sets would allow 
unambiguous determination of U. 

Having established a quantitative model for the SW 
dispersion allows to predict further quantities. We com- 
pute the non-interacting two-magnon dynamical struc- 
ture factor 5^^(k, w) probed by INS, and the two- 
magnon density of states (DOS) underlying, at (0, 0), 
Raman scattering and, at (tt, 0), K-edge RIXS jl4| . 
Although higher-order ma gno n interaction affect those 
two-magnon quantities |15l4l7| our results already allow 
several observations. Compared to IJHei, our predic- 
tions show the enhancement of a 500 meV peak in S"^^ 
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FIG. 4. (Color Online) (a-b) ^^^(k,^^) for IJHei with Jnn = 
^(tt,o)/'2Zc and the IbHub with the SCOC parameters, (c-d) 
energy lineshapes at the ZB points for IJHei (dashed red 
line) and IbHub (solid blue line), (e-f) Corresponding two- 
magnon DOS. (g-h) DOS profiles at (0,0) and (7r,0). (i-j) 
Single magnons wave-vectors contributing to the DOS peaks 
at (0, 0) and (tt, 0) in IJHei, (kd) in the SCOC IbHub. 



at (tt, 0) [Fig. 4(d)] which shows that attempts to ex- 
plain the reported INS lineshape at (tt, 0) from quan- 
tum effects must consider the full Hamiltonian presented 
here 0, [H, Also, along the ZB, the intensity of 

the one-magnon (transverse) excitations is constant so 
that the missing spin-wave amplitude observed by INS 
0, [H , 13, 21 1 does not result from further neighbor hop- 
ping. 

In the IJHei, the (0,0) two-magnon DOS peaks at 
4ZcJnn, corresponding to creating two spin- waves at the 
ZB. The peak in Raman Big spectra is found at 0.37 
eV for SCOC [13] corresponding to - 2.8Jnn- The re- 
duced energy was explained as due to magnon-magnon 
16|, but the peak- width could not be rc- 



mteractions 

produced. The large ZB dispersion that our model en- 
tails firstly imply that experiments should not be com- 
pared to a single Jnn, secondly it explains the Raman 
peak width as a range of energies from 2w(^^o) extending 



down to 2a; 



( 2 ' 2 ) ' 



where a maximum occurs because at 



this energy entire lines in one-magnon momentum space 
contribute [Fig. 4(k)]. Thirdly, it predicts a lower peak 
energy requiring weaker magnon interaction to match ex- 
periments. For a correct calculation of magnon interac- 
tions, we caution that in the cuprates, it is different one- 
magnon that contribute to the Raman peak [Fig. 4(k)], 



than in the IJHei [Fig. 4(i)]. The observation of a strong 
excitation at (tt, 0) in K-edge RIXS is also explained 
by our calculations, which demonstrate a concentration 
of DOS at exactly this wave-vector. Again the one- 
magnon states that contributes to this peak [Fig. 4(1)] are 
very different from the IJHei [Fig. 4(j)]. Thus, our re- 
sults reveal dramatic differences in the two-magnon con- 
tinuum, implying that INS, L3 or K-edge RIXS and Ra- 
man data must be interpreted using the full quantitative 
model derived here. 

Our model also provide insight into the electronic spec- 
tra, such as the bare-band dispersion used to extract 
the self-energy function from ARPES spectra. A self- 
consistent Kramers-Kroniger analysis of ARPES experi- 
ments on Bi2212 revealed t ~ 0.23 e V 1231 . However, with 
the discovery of a high-energy kink [2J[ at 0.4 eV in the 
nodal spectrum (also known as the waterfall feature), a 
similar analysis on optimally doped LSCO suggested 
that t ~ 0.48 eV. Our resuhs for LCO {t = 0.492 eV) 
support the second scenario. 

In summary we derived an effective spin Hamiltonian 
valid for any lattice and any hopping matrix element 
range. Using spin- wave theory with l/S'-corrections and 
three hopping t, t' and t" , we obtain accurate quantita- 
tive IbHub parameter sets for several parent compounds 
of the high-Tc cuprate superconductors. We predict or- 
dered moment, double occupancy, SW renormalization 
and 2-magnon spectra. From the non-interacting two- 
magnon S^^{q,u!) and DOS, we clearly demonstrate the 
necessity to include the extended exchange paths to inter- 
pret Raman and K-edge RIXS peaks. Furthermore, elec- 
tronic spectra such as ARPES could also be addressed 
using the same IbHub parameters. 
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